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ABSTRACT 

A generalization  of  the  Chandrasekhar  H-equation  is  solved  by  iteration. 

Such  equations  are  of  interest  in  heat  transfer. 

t 

SIGNIFICANCE  AND  EXPLANATION 

In  certain  mathematical  models  of  physical  phenomena,  there  may  be 
several  theoretical  solutions,  but  only  one  of  them  may  make  sense  in  the 
physical  context.  When  one  encounters  such  a model  which  cannot  be  solved 
in  closed  form,  the  question  naturally  arises  as  to  whether  some  numerical 
method  will  in  fact  produce  the  physically  significant  solution. 

In  certain  heat  transfer  problems  the  following  integral  equation  is 
encountered; 

(*)  H(x,w)  = 1 + H(x,a))  /“  J^(t)H(t,oj)dt,  x ^ 0 . 

In  (*),  (t)  is  the  Bessel  function  of  order  1,  w is  a parameter,  and 

h is  the  function  to  be  determined.  Of  the  several  possible  solutions  to 
(*) , only  one  is  of  physical  interest. 

In  this  paper,  we  show  that  (*)  may  be  solved  by  iteration  and  that  the 
solution  obtained  in  this  way  is  the  one  of  physical  inportance. 

AMS(MOS)  Subject  Classification:  80.45,  45G99,  45L05 
Key  Words;  H-equation,  iteration 
Work  Unit  No.  1 - Applied  Analysis 
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I . Introduction. 

Consider  the  integral  equation 

(1.1)  F(x,ui)  = 1 + (j  (t)  f(t,u))dt. 

0 ' 

In  (1.1),  is  a complex  parameter,  f is  the  function  to  be  found,  and  is  a measur- 

able function  on  (0,*). 

Under  certain  assumptions  on  it  is  knovm  that  (1.1)  has  a solution  H(x,u))  that  is 
analytic  in  w for  |u)|  <1.  This  solution  is  of  importance  in  many  applications.  Vte  call 

U the  physical  solution  to  (1.1).  The  question  considered  in  this  paper  is  the  following: 
When  can  H be  found  by  solving  (1.1)  by  iteration? 

/**  2 

t '^>{t)dt 

1 

<“  . In  this  Ccise  H is  continuous  for  |u|  ^1,  x ^ 1.  These  assumptions  are  natural  in 

/OO 

sC  (t)  dt  = — then 

1 

H may  be  found  for  |u|  < 1 by  an  iterative  method.  In  a more  general  setting,  Mullikin 
and  the  author  [7]  showed  that  H may  be  computed  for  lu|  ^ 1 in  this  way. 

In  many  cases  of  physical  importance,  however,  may  t)ecome  negative.  The  methods  of 
Bcwden  and  Zweifel  [1]  and  Rail  [9)  , then  give  the  best  known  results.  If,  however,  is 
not  integrable  on  (0,“),  these  methods  fail  as  the  integral  in  (1.1)  need  not  even  be  defined 
in  the  Lebesgue  sense . 

This  paper  is  motivated  by  work  of  Crosbie  and  Sawheny  [51,  [6].  They  consider  (1.1)  for 
(x)  = ^ Jj  (x)  , where  is  the  Bessel  function  of  order  one.  They  obtained  numerical 

results  that  indicated  that  H may  be  found  by  an  iterative  method  for  quite  general  . 

We  assume  •(  satisfies  the  following  conditions: 

(Al)  ; • Lp(0,0)  for  some  P , 1 ^ P < “ . 

(A2)  r 1^^^  |dt  <«  . 

0 

(A3)  Let  k(x)  “ lim  e t '^(t)dt;  then  k ^0,  and  k c Lj^  (-“,<»)  . 

N -*  » 0 

(A4)  / k(x)dx  ” 1 . 

••to 
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(A5)  lim  >p(t)dt  = 1/2  . 

N “ 0 

(A6)  If  f is  a measurable,  nonnegative  decreasing  function  on  (0,”)  and 

0 < sup  |x  f(x)|  < then  lim  i^(t)f(t)dt  > 0. 

0<x^l  N-*'"0 

•nie  reader  should  note  that  (A4)  is  merely  a normalization  condition  on  ^ . As  will 
become  clear,  (A5)  implies  (A4) . The  purpose  of  the  normalization  is  to  make  H analytic 
in  the  disc  |u|  < 1.  Our  assumptions  imply  that  the  integral  in  (1.1)  is  defined  for  non- 
negative decreasing  f by 

(1.2)  /"  ^(t)f(t)dt  = lim  /”  ¥>(t)f(t)dt  . 

0 N " 0 

For  e > 0,  let  C(e)  be  the  space  of  con^ilex-valued  functions  continuous  on  [e,") 
and  having  finite  limits  at  infinity.  C(e)  is  a Banach  space  under  the  sup  norm.  Moreover 
Dini's  Theorem  holds  for  C(e)  in  the  sense  that  if  a sequence  of  real  valued  functions  f^ 
in  C(e)  converge  monotonically  upward  to  f in  C(e)  pointwise,  and  lim  f^(x)  = limf(x) 

X -<■  " X-*  " 

for  all  n , then  f^  converges  to  f in  the  topology  of  C(e). 

We  will  show  in  section  II  that  H is  an  analytic  C(0)  valued  function  of  u for 
lu)|  <1  and  is  a continuous  C(e)  valued  function  of  u)  for  lull  ^1. 

Under  assun?)tions  (A1)-(A6)  we  prove  the  following  iteration  results. 


Theorem 

(2.1)  For  lu|  ^ 1 , 

X _>  0 

define 

Hp  (x,u)) 

= 1 

Hn^l(x.u3)  = 1 

+ loH  (x 
n 

fCO 

/ 

0 

H (t,u)(;' (t)dt 
x+t  n 

f n 

Then  for 

0 < 

e < 1 

(i) 

H 

n 

converges  to 

H ^ 

(7(0) 

uniformly  in 

U) 

for 

i<^i  1 

1 - e. 

(ii) 

H 

n 

converges  to 

H ^ 

(7(E) 

uniformly  in 

U) 

for 

i“l  1 

1 . 

Theorem 

(2.2) 

Let  0 < e 

< 1.  1 

u|  ^ 1 

, define  for 

X 

0 

Then 

Kg  (X,U) 

Vi'*’ 

= 1 

^00  f* 

“)  = ^ 
0 

'P(t)K  (t,V>)dtr 
n 

-1 

(i) 

K 

n 

converges  to 

H in 

C(o) 

uniformly  in 

OJ 

for 

|w|  < 

1 - e 

(ii) 

K 

n 

converges  to 

H in 

CU) 

uniformly  in 

(&> 

for 

W 1 

1 . 
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Implicit  in  the  statements  of  these  theorems  is  the  fact  that  all  integrals  exist  in 
the  sense  of  (1.2).  For  convenience  we  write,  for  g measurable, 

(1.3)  lim  g(t)dt  » / g(t)dt 

N + “ 0 * 

when  the  limit  on  the  left  exists. 


II  Iteration  Results 

Equation  (1.1)  is  intimately  connected  to  the  theory  of  Wiener-Hopf  equations  (8).  This 
connection  plays  a vital  role  in  the  proofs  of  our  results,  and  we  give  the  relevant  details 
here.  For  f e (-“,”),  f will  denote  the  Fourier  transform  of  f . 

Let  k(x)  be  as  in  (A3).  Assuit{>tions  (A3)  euid  (A4)  imply  that  for  |(i)|  < 1,  the  equation, 

/OO 

k(x-y) Y(y,(D)dy  = uk(x),  x > 0 , 

0 

has  a unique  solution  Y<x,<o)  e L^(0,<‘).  For  0 £ (i)  < 1,  Y ^ 0,  and  Y is  an  analytic 

Lj^  (0,") -valued  function  of  u for  |u)|  < 1.  We  write 

00 

(2.2)  = 51  u^Y  (x)  • 

n=l  ” 

In  (2.2),  Y (x)  >0  for  x > 0,  Y £ L (0,<“).  Then  the  unique  solution  to  (1.1)  which  is 
n — n i 

analytic  in  u for  |u|  <1  «md  C[0,<»]  valued  is 

(2.3)  H(x,o))  - 1 + y(ix,<a)  . 

For  |(d{  < 1,  we  have  the  factorization,  valid  for  real  X , 

(2.4)  (1  - u)k(X))  (1  + y(X,u))  (1  + Y(-X,ii)))  = 1 . 

From  (2.3)  and  (2.4)  we  conclude 

(2.5)  r — H(x,u)  < 0 for  x>0,  0<u<l, 

0 X — 

and 


(2.6) 


-7 — H(x,b))  > 0 for  x>0,  0<iD<l. 

<?  U)  — — 


(2.7) 

(2.8) 
(2.9) 


H(x,(ii)  = 1 + 5 (ix)  . |“|  < 1. 
n=l 

Y (ix)  > 0 for  n > 1 , x > 0 

n — — 

Y (i*)  e C[0,“]  for  n > 1 

n — 


Let  Z denote  the  class  of  nonnegative  decreasing  continuous  functions  f on  (0,”) 
which  satisfy 

(2.10)  sup  |xf  (x)  I < “ . 

0 < X < 1 
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For  t e Z define 

(2.11)  L f(x)  - ! -5-- ^(t)f  (t)dt 

' x+t 

* 

(2.12)  M f(x)  - i/>(t)f(t)dt. 

0 

T)ie  integral  defining  M is  a Lebesque  integral  for  each  fixed  x > 0 by  assua^tion  (Al). 

Assumptions  (Al)  and  (A6)  imply  that  Mf  and  Lf  are  in  C{0)  for  every  f e . 

We  note  that  for  0^u<l,  H is  in  We  may  rewrite  (1.1)  as 

(2.13)  H(u)  - 1 + u L(H(w))H(u)  . 

We  now  prove  the  main  lemma. 

Lenina  (2.1).  Let  c > 0,  then 

lim  _ H(u)  = H(l) 
u)  -♦  1 

exists  in  (7(e).  Moreover  H(l)  e Z , and  for  x > 0,  we  have 

(2.14)  H(l)  - 1 + L(H(1))H(1)  . 

Proof.  By  (2.4)  we  have  that 

H(0,u)  = (1  - = 1 + u / (t)H(t,u)dt  H(0,u). 

* 

Hence 

(2.15)  H'^(oa)  = (1  - u)^^^  + ID  M(H(o)))  . 

Consider  the  function  F given  for  x ^ 0 lay, 

(2.16)  F(x)  - / -^'(>(t)dt. 

Assuoptions  (A2)  and  (A6)  imply  that  F e CiO)  and  F > 0.  For  £ > 0,  x 0,  we  have 

F(x)  - F(*te)  - r (x^t)tx;'eW)  " ° • 

Therefore  F is  decreasing.  By  (2.13),  we  have 

(2.17)  h"^(x,u))  id  x F(x)  ^ id  x F(0)  . 

Hence,  for  x>e,0<iD<l,  H(x,id)  < — p”^  (0)  . Hence,  for  x > 0,  lim  H(x,id) 

— — — IDt  . 1 * 

ID  -*•  1 

« H(x,1)  exists  by  (2.6).  Moreover  H(x,l)  is  decreasing  in  x and  SV5>  |x  H(x,l)  | 

0<x  ^1 

< p”^(0)  < “.  Hence  M H(l)  e C(0)  by  assumption  (Al).  By  the  dominated  convergence  theorem 

and  assunption  (Al)  we  have,  for  x > 0 , 

M H(l)  (x)  - lim  _ M H(id)(x)  » h“^(x,id)  . 

ID  -►  1 

H(l)  e Z and  H(l)  satisfies  (2.14).  Dini's  theorem  then  implies  that  the  limit  in 
exists  in  Clt)  for  e > 0 . This  con^letes  the  proof. 


(2.18) 

Hence 

(2.18) 
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H(x,  |(»)|)  <_  H(x,l)  . 


Hence  for  e > 0 


Now  for  X > 0 , 


|h(x,u)) 


1 I l“rYn(ix) 


+ 1 


H is  an  cmalytlc  C(€)-valued  function  of  (d  for  lu|  < 1.  For  e > o,  H is  a continuous 
C(e)  valued  function  of  uj  for  |(a)|  < 1 . 

Vte  now  prove  the  main  results. 

Theorem  (2.1)  Assume  if  satisfies  (A1)-(A6).  For  |u|  1 , x ^0  define 


Hg(u)  = 1 

H , (uj)  = 1 + u)H  (iij)  L(H  (u) ) n > 0 . 
n+1  n n — 


Then,  for  0 < e < 1 , 

(i)  converges  to  H(u)  in  C(0)  uniformly  in  o)  for  |oj|  ^ 1 - e. 

(ii)  converges  to  in  C(e)  uniformly  in  u for  |u|  ^ 1 . 

Proof.  Let  0 < e < 1.  Osserve  that  if  g(x)  = 1 + J*  e"*^  G(t)dt,  and  G e L is  non- 

0 ^ 

negative , then  1 t Lg (x)  has  the  same  form  as 

(2.19)  Lg(x)  = /"e“’^[)c(t)  + /"  )c  (x+t)G  (s)ds]dt  . 

0 0 


The  definition  of  implies,  therefore  that  for  m ^1 

£ L,  (0,“)  , P > 0 , such  that  for  x > 0 , lool  < 1 , 

1 n ,ni  — — If  — 

2"-^ 

(2.20)  H (x,a))  = J u"*  /"  P (x)dx  . 


there  are  functions  P 

n,m 


It  is  easy  to  see  that 

(2.21)  P (x)  = Y (x)  , 1 < m < n . 

n,in  m — — 

vie  claim  that  P (x)  < y (x)  for  all  m and  n . We  prove  this  by  induction  on  n . For 
n,m  — m c-  j 

n = 0 the  result  is  clear.  Assume  (2.21)  for  n = N:  from  the  definition  of  H and  (2.19) 

n 

we  have, 

‘2-22)  Vl  m‘^’  = ^ P (t-s)  /”  )£(s+r)P  . (r)drds 

j,(l  > 1 

+ P„  „ 1 (t-s)lc(s)ds  + /“  lc(t+s)P„  , (s)ds 
Q N,m-1  ' N,m-1 

1 I Yj(t-s)  /“  )c(s+r)Y  (r)drdx 

j + (t=m-l  0 0 ^ 

j.i  1 1 

+ Y , (t-s))c(s)ds  + /”  )c(t+s)Y  , (s)ds 
0 0 

• Y (x)  . 
n 

The  last  equality  is  a consequence  of  the  fact  that 
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(2.23)  y (ix)  - [ Y.  (ix)L(  ^ , (ix)  + L(Y(i*))(x)  , 

in  • n ^ D ^ in~i 


j.i  i 1 

which  is  an  easy  consequence  of  (1.1). 
A similar  argument  shows  that 


<2-24)  „ 1 m • 

n,in  — n+1  »in 

Hence  we  have,  for  n,k2.0»  x^O,  |u(^l, 

OD 

(2.25)  |h  (x.u)  - H (x,u))l  < H(X,  lul)  - H (x,|ul)  < I Y„(ix)|u|'”  . 

" m-n+1  “ 

For  X ^ e , |(i)l  1 we  have 

® n 

(2.26)  I Y (ix)  liij"  1 H(x,l)  - I y 

mi«n+l  m-1  m 

The  right  side  of  (2.26)  goes  to  zero  as  n becomes  large  uniformly  for  x >^  £ , |ii)|  ^ 1- 

This  proves  (ii). 

For  x^O,  |u|^l-e  we  have 

^ n 

(2.27)  I Y(ix)  Iwl”  H(x,l-e)  - I Y(ix)|l-e|”  . 

m-l 

m-n+1 

The  right  side  of  (2.27)  converges  to  zero  as  n becomes  large  uniformly  for  x ^ 0 , 
|u|  ^ 1-e  . This  proves  (i)  and  con^letes  the  proof. 

Theorem  (2.2)  Let  0 < e < 1,  define,  for  juj  ^1,  x 0 , 

Kq  (x , u)  « 1 

K , (x,u)  - [1  -u  L(K  (*,u))(x)l  ^ , n > 0 . 
n+1  n “ 


(i)  K converges  to  H in  C(C)  uniformly  in  oj  for  |aj|  < 1-e  . 

n 

(ii)  K converges  to  H in  C(e)  uniformly  in  u for  |u|  < 1 

n 

Proof.  Note  that  Hjj(x,u)  « Kp(x,ai)  ^ H(x,u)  . An  induction  argument  similar  to  that  used 

in  the  proof  of  Theorem  (2.1)  will  show  that,  for  x>0,  Of_ii)f_l. 

(2.28)  K (x,a))  < H(x,a))  , 

n — 


(2.29)  K (x,ii))  ll-u)  L(H  .(•,id)(x))  >H  (x,a))  . 

n n-1  ~ n 

As  in  the  proof  of  Theorem  (2.1)  we  also  have  for  lu)|  f_l,  x > 0 , n ^0,  m^O, 

(2.30)  (K^^^(x,a))  - K^(x,u))|  ^ (x , | oi | ) - K^(x,|u)(). 

This  completes  the  proof. 

Systems  of  nonlinear  equations  similar  to  (1.1)  are  of  interest  in  the  kinetic  theory  of 
gases  [2],  [4].  It  is  possible  that  the  methods  of  this  paper  will  generalize  to  that  setting. 
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